Frustration-induced r] inversion in the S=l/2 bond-alternating spin chain 
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We study the frustration-induced enhancement of the incommensurate correlation for a bond- 
alternating quantum spin chain in a magnetic field, which is associated with a quasi-one-dimensional 
organic compound F5PNN. We investigate the temperature dependence of the staggered suscepti- 
bilities by using the density matrix renormalization group, and then find that the incommensurate 
correlation becomes dominant in a certain range of the magnetic field. We also discuss the mech- 
anism of this enhancement on the basis of the mapping to the effective S=l/2 XXZ chain and a 
possibility of the field-induced incommensurate long range order. 
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The one-dimensional (ID) S = 1/2 antiferromagnetic 
bond-alternating spin chain has been an important is- 
sue in the condensed matter physics, since it exhibits 
some typical quantum effects. The magnetic susceptibil- 
ity clearly reflects the existence of the dimer spin gap. A 
more important aspect is that the Tomonaga-Luttinger 
(TL) liquid, which is an undoubtedly essential concept 
in ID quantum critical systems, is realized between the 
vanishing field of the spin gap Hd and the saturation 
field Hc2, where the low-energy behavior of the system 
is characterized by the TL exponents rj^ {riz) associated 
with the correlation function for the transverse (longitu- 
dinal) staggered mode (see Eq. l|2Il) 

A good example of such a bond-alternating spin chain 
is a quasi-lD organic compound F5PNN |^. However, a 
recent precise analysis of F5PNN suggests that the frus- 
tration effect, which is recently attracting considerable 
attention, induces various anomalous properties (J, 0, Q ■ 
In the magnetization process, the bond alternation ra- 
tio a (defined in Eq. (^) shows a crossover from 0.4 in 
the low field region to 0.5 in the high field region j^f. 
Moreover, the temperature dependence of the NMR re- 
laxation rate in a magnetic field exhibits anomalous en- 
hancement of the TL exponent ry^; in contrast to the 
usual behavior of the TL exponents r]x < rjz for the sim- 
ple bond-alternating system 01 j > ??z may be realized 
in a certain range of the magnetic field 0, Q , which we 
shall call "ry-inversion" in this paper. 

Motivated by such interesting experimental sugges- 
tions of the frustration effect, we study the bond- 
alternating spin chain with the frustration in a magnetic 
field H: 

H = J ■ <5'2i-|-l + 0!S2i+l ■ 821+2] 

i 

+ J'Y,S,-S,+2-HY,St, (1) 

i i 

where S is the S = 1/2 spin operator and J' is the frus- 



trating coupling. In the following we use the normal- 
ization J = 1 for simplicity. Particular importance of 
this model is that it enable us to capture the interesting 
physics cooperatively generated by the magnetic field and 
frustration. In fact the frustration-induced plateau for- 
mation has been studied intensively j2i Isl , and recently 
a remarkable enhancement of the incommensurate cor- 
relation rjx > rjz has been suggested by the numerical 
diagonalization analysis in some intermediate magnetic 
field 9]. However, a systematic investigation of the frus- 
tration effect for the TL liquid behavior in the magneti- 
zation curve is essentially difficult and a detailed study 
is clearly desired for more thorough understanding of its 
impact on the TL exponents. In addition, a precise anal- 
ysis for the 7y-inversion provides an essential view point 
for the material physics; as an interesting consequence of 
such enhancement of the incommensurate correlation, a 
novel type of incommensurate order can be induced in 
the magnetic field through the inter-chain coupling. 

In this paper, we reveal the effect of the frustrating 
coupling J' on the observable quantities, using the fi- 
nite temperature density matrix renormalization group 
(DMRG) d m We first calculate the magne- 

tization curve and further investigate the temperature 
dependence of the staggered susceptibility xa- (X||) Per- 
pendicular (parallel) to the uniform magnetic field. We 
then find that the ry-inversion actually occurs in a certain 
range of the magnetic field, where the gap formation at 
the M = 1/4 plateau plays a crucial role. Here M is the 
magnetization per one spin. Finally we discuss the possi- 
bility of the field induced incommensurate order assisted 
by the frustration. On the basis of the obtained phase 
diagram, we make a comment on the NMR experiment 

of F5PNN. 

In order to analyze the crossover of a in F5PNN, we 
first calculate the magnetization process by the finite 
temperature DMRG with the retained bases number 64. 
Figure n shows the comparison of the obtained curves at 
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T = 0.085. The parameters in the figure correspond to 
the F5PNN experiment in Ref. Q. We can clearly see 
that the curve of a = 0.45 with the frustrating coupling 
J' — 0.05 explains the crossover of the magnetization 
process from a ~ 0.4 in the low field region to a ~ 0.5 in 
the high field region. This gives a clear evidence of the 
frustration effect in F5PNN. 




FIG. 1: Magnetization curves obtained with the DMRG. 

Let us next discuss the TL exponents 77^ and rjz, the 
precise definitions of which are given by power-law decay 
of the spin-spin correlation functions at zero tempera- 
ture: 

{S^S^^)~M'^ - cos{2kFr)r-'^% (2) 

where kp = 7r(l/2 — M). It is well known that these 
exponents provide the essential information for the ob- 
servable quantities. For instance, the low temperature 
behaviors of the staggered susceptibilities are character- 
ized by these exponents [l^ : 

Xi(r)^T-(2-''^) ;t||(T)^r-(2-')^). (3) 

Since the relation rjx'qz = 1 is satisfied for the TL liquid, 
the smaller exponent yields the dominant spin fluctuation 
in low temperatures. Since the bond-alternating chain 
usually has rjx < rj^, x± shows the stronger divergence 
in the T — > limit. However, if the incommensurate 
spin correlation along the magnetic field is enhanced, X|| 
can be dominant in the low temperature region. In or- 
der to see how the frustrating coupling J' affects the in- 
commensurate correlation, we directly calculate the stag- 
gered susceptibilities y i (T) and X||(^) with the finite 
temperature DMRG [liT llll ll^. which is the most reli- 
able method to investi gate an infinitely long chain with 
the frustration [l^ Hslllij. In actual computations, we 
selectively employ the infinite size and finite size algo- 
rithms of the DMRG. For a calculation of x±, we per- 
form the infinite size DMRG with a weak commensurate 
staggered field along the x direction; x± is obtained as a 
numerical derivation of the staggered field. In contrast, it 
is generally difficult to treat directly the 2fci?-oscillating 



magnetic field with the finite temperature DMRG. Thus 
we start with the linear-response formula: 

X||(T)= J2 r dr{S^-{T,r)S^ {0,0)), (4) 

where /3 is the inverse temperature, Tc is a cut-off for 
the real space direction, and (r, r) is the spin operator 
in the Heisenberg representation at an imaginary time r 
and position r. By using the finite size algorithm of the 
DMRG for the quantum transfer matrix T, we calculate 
the correlation function (iS^(t, r)5^ (0, 0)). After obtain- 
ing the maximum eigenvalue of T and the corresponding 
eigenvector li/^max), we calculate {S^ {T,r)S^ (0,0)) as 

{S (r, r)S (0, 0)) = . (5) 

Here 5*"^. is the spin operator at a position (r, r) on the 
checkerboard lattice obtained via the Suzuki- Trotter de- 
composition [T^ . In the numerical computation, the in- 
tegral for T in Eq. Q is replaced by the summation with 
a finite imaginary-time-slice e = P/N, where N is the 
Trotter number. In the following results, we have set 
Tc = 200 and iV = 80, and confirmed the sufficient con- 
vergence with respect to Vc and N. Here, it should be 
noted that, especially at M = 1/4, the consistency of the 
finite size algorithm can be directly checked with the nu- 
merical derivative of the infinite-size DMRG. This is be- 
cause the renormalization process for the quantum trans- 
fer matrix is compatible to the periodicity of kp = n/2 at 
M = 1/4. We have also confirmed that the both results 
are in good agreement with each other. 

In Fig. we show comparisons of x± ^-^d X|| 
for {a,J',H) = (0.45,0.05,1.16) and {a,J',H) = 
(0.45,0.15, 1.19). In both cases, the corresponding mag- 
netizations at the zero temperature are slightly lower 
than M = 1/4. We can see that xj_ for J' = 0.05 is 
always larger than xy, implying that the commensurate 
fluctuation is still dominant. The estimated exponents 
are rix = 0.78 and 77^ = 1.3, which yields ri^riz = 1.0. In 
contrast, for J' = 0.15, the susceptibilities clearly show 
the crossover around T 0.1; as the temperature is de- 
creased below T ~ 0.1, the divergence of x± is reduced, 
while that of x\\ is enhanced. The TL exponents in the 
low temperature region are extracted as r/a, = 1.3 and 
r]z — 0.8, which also gives rix'i]z 1.0. We can thus verify 
that the ry-inversion is actually realized by the frustration 
effect. In addition, we have found that, as T — > 0, x\\ iii 
the M = 1/4 plateau diverges exponentially while x± 
converges to a finite value. This fact suggests that the 
appearance of the plateau gap plays an important role 
for the 77 inversion. 

In order to discuss the microscopic origin of the 77 in- 
version, let us recall that the Hamiltonian (^) around 
the plateau can be mapped to the effective S=l/2 XXZ 
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FIG. 2: Staggered susceptibilities X|| X± {^) for 

{a,J',H) = (0.45,0.05,1.16) and (b) for {a, f , H) = 

(0.45,0.15,1.19). The dotted (thin solid) lines show the re- 
sults of power- law fitting for X-l(X||)- 



chain where the anisotropy and magnetic field of 



the effective model 



by A 



1 2J'+a 



and 



are given ^ - 2\2J^~ 
-ffoff — H — 1 — (2 J' + q;)/4 respectively. The exact crit- 
ical exponents of the XXZ chain in a magnetic field can 
be derived from the Bethe ansatz integral equation for 
the dressed charge Since the rjz of the S=l/2 XXZ 
chain is a monotonous increasing function of |fleff|, it is 
sufficient to investigate rjz at iJeff = for the purpose 
of understanding the appearance of the ry-inversion. For 
the XXZ model, the following fact is well known; as far 
as A < 1, the XXZ chain is in the critical TL phase with 
r]z > 1, while, for A > 1, the excitation gap is opened 
and, at the same time, rj^ < 1 appears only near the 
zero magnetic field. We therefore find that the criterion 
for rix > rjz is equivalent to the one for the gap forma- 
tion. Turning to the original model, we can see that the 
condition for rjx > rjz is deduced as a/6 < J' < 3a/2, 
which is satisfied by the parameters used in Fig. |21 (b). 
Since the effective anisotropy A becomes large as J' ap- 
proaches a/2, we can also understand that the region 
Vx > rjz extends as J' — > a/2. Although the mapping 
to the XXZ chain is based on the perturbation theory 
from the a, J' — > limit, the precise study by the level 
spectroscopy method ^ gives a good estimation for the 
1/4 plateau, where the effective XXZ model picture is 
basically maintained In fact, rjz of the effective XXZ 



model for {a,J',H) = (0.45,0.15,1.19) is calculated as 
rjz = 0.735, which is consistent with the fitting result in 
Fig- El (b). Here we should note that, for a = 0.45, the 
1/4 plateau at T = emerges for 0.08 < J' < 0.5 0, im- 
plying that F5PNN is located at a subtle position near 
the plateau phase boundary. 

On the basis of the results mentioned above, let us dis- 
cuss the effect of the inter-chain interaction Jjnt, which 
induces the 3D long-range order corresponding to the 
dominant spin correlation [20j . Indeed the sharp peak of 
the specific heat was observed for F5PNN in a magnetic 
field For the weakly frustrating bond-alternating 

chain, the antiferromagnetic long-range order is realized 
for Hci < H < Hc2- However, when the strong frus- 
trating coupling J' induces the r/-inversion, the 3D order 
should also change from the ordinary Neel order perpen- 
dicular to H into the incommensurate order along H at 
some field. In order to clarify whether such a change of 
the order occurs or not, we calculate the theoretical phase 
diagram of the frustrated bond-alternating chain within 
the inter-chain mean field approximation combined with 
the DMRG HI HI; using x\\{T) and x±iT) obtained by 
the finite temperature DMRG, we can extract the phase 
boundary from the equation x-yi'I'c) — {zJint}^^ [23j |. 
where z is a coordination number and 7 G_L or |j. Here 
we assume that the inter-ladder coupling is weak and is 
not frustrating. 

In Fig. 13 we show the H-T phase diagrams for a — 
0.45 with J' = 0.05 corresponding to F5PNN, and with 
J' = 0.15. In this figure, the phase boundary for a fixed 
zJint is indicated by a line. For a = 0.45 and J' = 0.05, 
the transition to the commensurate order always appears, 
reflecting the fact x± > X\\ in Fig-2 (a). Moreover, we can 
see that the phase boundary for F5PNN is well re- 
produced by the theoretical curve for zJint = 1/20. The 
estimated coupling J/ fcs = 6K is also consistent with the 
one obtained in Ref. This fact provides another veri- 
fication for J' = 0.05 determined from the magnetization 
curve. If we assume z = 4, the inter-chain coupling is es- 
timated as Jint/^B — 0.07K. As J' is increased, the usual 
commensurate fluctuation is suppressed around 7? ~ 1 
and, at the same time, the fluctuation associated with 
rjz is increased. Accordingly the phase boundary of the 
commensurate Neel order is shifted to the low temper- 
ature side and, as far as x\\ — ^/zJhit is satisfied, the 
incommensurate order can appear in the left side of the 
bold line in Fig. 3 (b), which becomes identical to the 
commensurate-incommensurate transition line for suffi- 
ciently large Jjnt. As was seen before, the 77 inversion si- 
multaneously appears with the plateau formation. Thus 
we can understand that the incommensurate order devel- 
ops in the intermediate field region. 

Finally, we want to make a comment on the NMR 
relaxation rate of F5PNN, for which the theoretical H- 
T phase diagram does not show the ry-inversion. Since 
the temperature region used for the estimation of rjz in 
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Ref. |5|, Ig is close to the transition temperature, the ex- 
perimentally observed crossover oirjz for F5PNN is possi- 
bly due to the anomaly originating from the transition to 
the commensurate Neel order. However the 77 inversion 
certainly occurs in the large J' region. 
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FIG. 3: Phase diagram determined by the inter-chain mean 
field approximation (a) for (a, J') = (0.45, 0.05) and (b) for 
{a, J') = (0.45,0.15). Pluses show the critical temperature 
tJ{H) of F5PNN with J/ks = 6K. The bold line in (b) 
indicates the X|| = X-L line. Dotted lines represent the critical 
fields {Hci and i/c2), and the lower (upper) edge of the plateau 
Hpi {Hp2) for the single chain. 

To summarize, we have studied the bond-alternating 
spin chain with the frustrating interaction in a magnetic 
field. What we want to emphasize in the present study is 
that the combination of the magnetic field and frustration 
gives rise to various exotic phenomena. We have actually 
shown that the frustration explains the crossover of a in 
the magnetization curve of F5PNN. We have also clari- 
fied that the ry-inversion of the TL exponents is induced 
in a certain range of the magnetic field by the frustra- 
tion. The mechanism of the 77-inversion is explained on 
the basis of the effective XXZ model associated with the 
1/4 plateau formation. Moreover, we discussed the pos- 
sibility of the field induced incommensurate order. Al- 
though the nature of this transition between commensu- 
rate and incommensurate orders is not clear within the 
mean field theory, we believe that we could elucidate the 
importance of the incommensurate correlation. Recently, 
several frustrated spin systems have been studied, where 



a low-e nerg y effective XXZ model often works success- 
fully 0, 0, ■ I'^ particular the similar enhancement 
of r]x is also reported in Ref. [13 • This suggests that 
the mechanism for the 77-inversion can be applicable to a 
class of realistic frustrating spin systems. Then we hope 
that the present theory stimulates further researches on 
such interesting physics caused by the incommensurate 
correlation. 
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